Passive scalar intermittency in compressible flow 
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A compressible generalization of the Kraichnan model (Phys. Rev. Lett. 72, 1016 (1994)) of passive 
scalar advection is considered. The dynamical role of compressibility on the intermittency of the 
scalar statistics is investigated for the direct cascade regime. Simple physical arguments suggest that 
an enhanced intermittency should appear for increasing compressibility, due to the slowing down 
of Lagrangian trajectory separations. This is confirmed by a numerical study of the dependence of 
intermittency exponents on the degree of compressibility, by a Lagrangian method for calculating 
simultaneous A'^-point tracer correlations. 
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In the last few years, much effort has been devoted to 
the study of statistical properties of scalar quantities ad- 
vected by random flows with short memory. Remarkable 
progress in understanding intermittency and anomalous 
scaling has been achieved ||l|-Q for the Kraichnan model 
PI of passive scalar advection by random, Gaussian, in- 
compressible and white-in-time velocity fields. A crucial 
property of the model is that equal-time correlation func- 
tions obey closed equation of motion. Analytical treat- 
ments are thus feasible, and the identification of a gen- 
eral mechanism for intermittency has been established. 
Its source has been found in zero modes of the operators 
governing the Eulerian dynamics of A^-point correlation 
functions ||,|§. 

Concerning numerical studies of the Kraichnan model, ef- 
ficient Lagrangian methods have been recently proposed 
P,p| and thanks to them both the limits of the vanish- 
ing of intermittency corrections, for which perturbative 
predictions are available M-M, and the non-perturbative 
region, have been successfully investigated |||,g] . 

A compressible generalization of the Kraichnan model 
has been recently proposed 110-02] and the existence 
of very different behaviors for the Lagrangian trajecto- 
ries, depending on the degree of compressibility, has been 
shown analytically fl^Jllll . For weak compressibility, the 
well-known direct cascade of the passive scalar energy 
takes place. This is associated, from a Lagrangian point 
of view, to the explosive separation of initially close tra- 
jectories [pyi3|], a feature characterizing the direct en- 
ergy cascade for the incompressible Kraichnan model as 
well. On the contrary, when the compressibility is strong 
enough, particles collapse: both non intermittent inverse 
cascade of tracer energy exciting large scales and suppres- 
sion of the short-scale dissipation occur [Q . The relation 
between intermittency and compressibility is the main is- 
sue of the present short communication. 
As already highlighted [^ 11|, because compressibility 
inhibits the separation between Lagrangian trajectories, 
the resulting scalar transport slows down and scaling 
properties may be affected. Our remark here is that 
the slowing down of Lagrangian separations plays an es- 



sential role in characterizing intermittency in the direct 
cascade regime. This can be easily grasped from the fol- 
lowing considerations. In the direct cascade regime, typ- 
ical trajectories are stretched, whereas contractions are 
rare and thus affect only the extreme tails of the pdf 
of scalar differences. Furthermore, within a Lagrangian 
framework, scalar correlations are essentially governed 
by the time spent by particles with their mutual dis- 
tances smaller than the integral scale of the problem. 
The stretching process, typical of the direct energy cas- 
cade, is thus intermittent because contracted trajectories 
cause strong fluctuations of the time needed to reach the 
integral scale. When compressibility is present, even if 
weakly, trapping effects are amplifled due to the slowing 
down of Lagrangian separations. It then follows that the 
dynamical role of collapsing trajectories increases for in- 
creasing compressibility, and the same should happen for 
the intermittency. It is worth noting that the trapping 
mechanism, enhanced by the compressibility, works in 
the same direction as that induced by lowering the spa- 
tial dimension d: it is indeed observed perturbatively S] 
that when d is reduced an increased intermittency arises, 
a fact corroborated by numerical evidences f^ compar- 
ing results of the incompressible Kraichnan model in two 
and three dimensions. These considerations will be here 
quantitatively supported by numerical simulations. 

The compressible generalization of the Kraichnan 
model is governed by the equation (for the Eulerian dy- 
namics) 

dtOir, t) + t;(r-, t) • V 0(r, t) = KV''e{r, t) + /(r, t), (1) 

where, as for the incompressible case, the velocity and 
the forcing are zero mean, Gaussian independent pro- 
cesses, both homogeneous, isotropic and white-in-time. 
The velocity is self-similar, with the 2-point correlation 
function: 

{v^{r,t)vp{r',t'))^5{t-t') [d% - d^p{r - r')] , (2) 

where dai3{r), the so-called eddy-diffusivity, is fixed by 
isotropy and scaling behavior along the scales: 



dap{r) 



{[A+{d + i-\)B]5^p^i[A~B]'^] 



(3) 



where d is the dimension of the space. 
The degree of compressibility is controlled by the ratio 
p = C'^/S'^, being S^ = A+ [d- 1)B oc ((Vw)^) and C^ ee 
A ex ((V • v)"^), which satisfies the inequality < p < 1. 
The statistics of the forcing term is defined by the 2-point 
correlation function 



{f{r,t)fir',t'))=Sit-t')x{\r-r'\), 



(4) 



where x is chosen nearly constant for distance \r — r'| 
smaller than the integral scale L and rapidly decreasing 
for r ^ L. 

It is worth remarking that equation (m) physically de- 
scribes the evolution of a tracer, that is a quantity which 
is conserved along the Lagrangian trajectories in absence 
of diffusivity and forcing. To characterize the advection 
of a density, one should consider the equation 

dtpir, <) + V • [vir, t)p{r, t)) = ^V'p{r, t) + /(r, i), (5) 

which in the ideal case (k = 0, / = 0) enjoys the conser- 
vation of the total mass. The density advection equation 
has also a wide realm of physical applications and should 
deserve a detailed study in its own, as well as a specific 
numerical approach. Hereafter we shall limit ourselves to 
the case of tracer advection ruled by (|^). 

Exploiting the i5-correlation in time, equations for the 
even scalar correlations (odd correlations being triv- 
ially zero) in the stationary state, can be deduced 
|l4| ; for the generic A^-point correlation function C% = 
{9{ri) ■ ■ ■ 0{rN)) the expression reads: 



i<j 



... .0{rN)) , 



(6) 



with rij = ri — rj, and AIat is the differential operator 
given by: 



Mn = 

y^ dal3{rn-rm)Vr 

l<n<m<N 



l<n<N 



(7) 



As for the incompressible case, this model has a Gaus- 
sian limit for ^ — > 0, and the perturbative expansion 
at small ^'s can be done as in Ref. |g]. Accordingly, 
the calculation performed in the weakly compressible 
case (i.e. p < d/^^) corresponding to the direct cas- 
cade regime leads (see Ref. |0]) to the expression for 
the intermittent correction A^, to the normal scaling 
exponent (2 — ^)N/2 of the N-point structure function 
S%{r) = {[d{r) - e{0)]^) oc r(2-«)^/2-A?,. namely: 



A 



N 



N{N-2){l + 2p) 
2{d-2) 



i+oie)- 



(8) 



The perturbative approach gives thus a first clue that 
compressibility works to enhance intermittent correc- 
tions. We are however interested in checking that this 
is a general and robust feature associated to compress- 
ibility and thus that it is present for generic ^. This 
problem is not accessible by perturbative techniques; nu- 
merical methods are generally needed to investigate it. 
With this purpose in mind, we have developed a new La- 
grangian numerical method (a different viewpoint with 
respect to the one in Ref. g), where the strategy is now 
formulated in terms of a first exit time problem pq l . 

The method consists in the Montecarlo simulation of 
Lagrangian trajectories according to the stochastic dif- 
ferential equation 



(9) 



Virn,t) + V2K.Wn , 



where the w„ are independent Wiener processes. The 
evolution of the probability PN{t,x\to,Xo) that the N 
Lagrangian tracers have a configuration x = (ri , . . . , tm) 
at time t given their initial configuration xq at time ^o is 
ruled by the Fokker-Planck equation 



d_ 
di 



PN{t,x\to,Xo) + MNix)PN{t,x\to,Xo) =0, (10) 



where the operator A^^ is the adjoint of (H). As a con- 
sequence of ( [lO|) the probability obeys also the backward 
Kolmogorov equation 

d 

—-PN{t,x\to,Xo) + MN{xo)PN{t,x\to,Xo) = . (11) 
OTo 

We now introduce the Green function 

/•OO 

G{x,xo)= dt PN{t,x\ta,xo) , 

Jto 

which enjoys the following properties 



(12) 



M*j^{x)G{x,Xq) = ~S{X - Xq) 
Mn{Xo)G{x,Xo) = ~S{X - Xq) 



(13) 
(14) 



Let us define the characteristic size of a configuration 
of N particles as R{x) = [{J2,<j \r^ - r^f )/(iV(iV - 
l)/2)]^/^ . We now impose Dirichlet (absorbing) bound- 
ary conditions at R{x) = L ^ R{xo), and compute nu- 
merically the first exit time from the volume of configu- 
ration space limited by the boundary, which is expressed 
in terms of the Green function as (see e.g. [Em) 



Tl{xo) 



dx G{x,Xo). 



(15) 



lRix)<L 

A trivial consequence of the property (O) is that 

Mn{xo)Tl{xo) = -1, (16) 

an equation whose structure resembles that of (pi); indeed 
we can conclude, similarly to what happens for correla- 
tion functions (e.g. [|||]), that Tl{xq) must amount to 



the sum of an inhomogeneous solution plus a linear com- 
bination of zero modes fj of the operator Mn' 

nixa) = J2 C.L^-"^ fj{xQ) + inhomog. term , (17) 



where the explicit dependence on L has been extracted 
taking advantage of the scaling properties of A^at, aj 
is the scaling exponent of the zero mode fj and Cj is 
a constant independent of L. Among the non trivial 
zero-modes fj, only the functions which depend on all 
the coordinates can contribute to the A^-th order struc- 
ture function. We would like to extract this contribution 
leaving aside all the others: it is easy to realize that this 
result can be achieved performing a linear combination of 
the exit times with different initial conditions. This op- 
eration will remove also the inhomogeneous term. If we 
denote with Vi(p) the operator acting on the functions of 
A^-particles coordinates as Vi(p)F(ri, . . . , r^, . . . , rjv) = 
F(ri, . . . , r, -t- p, . . . , rA,) - F(ri, . . . , r,, . . . , rjy) we will 
have 



^n{L)^\1^,{p)Tl{x^)^L 



7-Ca 



(18) 



where Cat = (2 — £,)N/2 — A^ is the scaling exponent of 
the structure function S^{r) ~ r'»™. Whenever Xq — 0, 
due to the simmetry of the fj 's under exchanges of parti- 
cles coordinates, the expression for I]jv(i) takes a simple 
form, which, for example, for A^ = 4 reads as "S^lL) — 
2 mo, 0, 0, 0) -8Tl{p, 0, 0, 0) + 6Tl{p, p, 0, 0). 

Summarizing: the numerical method consists in the 
Monte-Carlo simulation of Lagrangian trajectories of A^ 
particles advected by a rapidly changing velocity field, ac- 
cording to the Fokker-Planck equation (nw ; average first 
exit times outside a volume of size L are computed for 
different arrangements of the initial conditions, and then 
linearly combined according to dlq) in order to extract 
the scaling exponent Cat. 

As a final remark, the numerical method here em- 
ployed can be viewed as a merging of the two Lagrangian 
methods introduced by Frisch, Mazzino and Vergassola 
in Ref. S and by Gat and Zeitak in Ref. M. Namely, 
it borrows form the first one the idea of subtracting exit 
times of different initial conditions to extract the only 
zero mode that contributes to the structure functions, 
while inherits from the second the spirit of working with 
particle configurations (shapes). The advantages of the 
present method with respect to ||^ mainly reside in the 
evaluation of first exit times rather than of residence 
times, a fact which substantially reduces the computa- 
tional cost. 

We present the numerical results obtained for the scal- 
ing of the fourth-order structure function S4^{r;L) = 
{{9{r) — 9(0))'^) in three dimensions. As previously men- 
tioned, when the dimension d of the space is lowered 
fluctuations increase and as a consequence the number of 
realizations needed to have a clean scaling grows as well; 



the addition of compressibility further enhances this ef- 
fect. For the first numerical experiments with the new 
method, we have thus opted for d = 3. 
The method has been tested performing the analysis of 
the incompressible limit p — Q for different values of ^: 
the anomaly A4 — 2C,2 — C4 has always been found to 
be compatible with the results presented in Refs. [^,^. 
The computation of S2 {L) - which can be evaluated ana- 
litically - has provided another stringent test for this 
method. 

Varying the degree of compressibility p, we have studied 
in the direct cascade regime the connection between the 
slowing down of Lagrangian trajectories and intermit- 
tency at the two distinct values ^ = 0.75 and ^ = 1.1. 
Notice that for these two values of ^, the condition 
(p < d/^^) for the direct cascade of energy to take place 
|ll| is verified for the entire range of values < p < 1 
of the compressibility. Different motivations account for 
this choice; first of all we avoided the region of ^ close to 
(7 — > 2) where capturing the subdominant anomalous 
exponents is numerically expensive, and furthermore the 
results are known from perturbative expansion. Second, 
when ^ is close to 2 (7 — > 0) non local effects are very 
strong and the range of values of p (i.e. p < d/^"^) per- 
taining to the direct cascade is narrower. 

In Figs. and g are shown the behavior of S4(i) for 
the two values of S, under consideration and for different 
values of p, which all display a fairly good power law scal- 
ing. According to the relation (n8) the scaling exponent 



is 7 - C4 



-7 -|- A4, so that the curves become flat- 



ter and flatter as the anomaly grows . It is thus evident 
from our results that when compressibility increases, the 
intermittent correction to the normal scaling grows as 
well. 
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FIG. 1. A log-log plot of Y,a{L) for ? = 0.75. (a): p = 0; 
(b): p = 0.25; (c): p = 0.5; (d): p = 0.75. Separation 
p = 2.7 X 10~ , diffusivity k = 2.3 x 10^"", number of realiza- 
tions ranging from 20 x 10® (case (a)) to 30 x lO" (case (d)). 
Solid lines represent the best fit power laws. 
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FIG. 2. As in Fig. 1, for C 
K = 2.5x 10"^ 



1.1 and diffusivity 



lationship between intermittency of passive scalar statis- 
tics and compressibility of the advecting velocity field. 
This result can be easily understood from the Lagrangian 
viewpoint. Intermittency arises whenever the particles 
experience long periods of inhibited separation: since 
compressible flows are characterized by the presence of 
trapping regions, an enhancement of intermittency can 
be reasonably expected. The validity of this argument 
has been assessed by means of a numerical Lagrangian 
method. 

We acknowledge innumerable discussions on the sub- 
ject matter with M. Vergassola. Simulations were per- 
formed in the framework of the SIVAM project of the 
Observatoire de la Cote d'Azur. Part of them were per- 
formed using the computing facilities of CINECA. 



Notice that ratio between S4 and the dominant con- 
tribution to each term of the sum scales as L"^''. As a 
consequence, small values of ^ (which correspond to large 
values of C4) require a larger amount of statistics to make 
the subdominant contribution emerge. This is the reason 
for which the scaling region for ^ = 0.75 is smaller than 
that for ^ = 1.1. 

Finally, our results are summarized in Fig. which 
shows the anomaly 2(^2 — Ci ^s the compressibility factor 
p for £_ — 0.75 (squares joined by a dot-dashed hue) and 
^ = 1.1 (circles joined by a dashed line). As in Ref. [||, 
the error bars are obtained by analyzing the fluctuations 
of local scaling exponents over octave ratios of values for 
L, a method which gives a very conservative estimate of 
the errors. The effectiveness of the first exit time compu- 
tation is somehow balanced by the need of a huge number 
of realizations to achieve a satisfactory statistical conver- 
gence. This drawback is particularly visible for large L, 
where the signal is rather noisy. 



1^ 0.50 



0.00 0.25 0.50 0.75 

Compressibility factor 

FIG. 3. The anomaly 2^2 —C4 for the fourth-order structure 
function, for ^ = 0.75 (squares joined by a dot-dashed line) 
and ^ = 1.1 (circles joined by a dashed line). 
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In conclusion, we have shown in the context of the 
Kraichnan compressible model that there is a tight re- 



